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Abstract 



We study the initial value problem of the thermal-diffusive combustion system: 
^ , ui^t = ui,x,x—ui'u\, U2,t = du2,xx+ui'>J'2, X £ , for non-negative spatially decaying 

I initial data of arbitrary size and for any positive constant d. We show that if the 

initial data decays to zero sufficiently fast at infinity, then the solution (ni,M2) 
converges to a self-similar solution of the reduced system: ui^t = ui^xx—uiU2,U2,t = 
du2,xxi in the large time limit. In particular, ui decays to zero like 0{t~2-^), where 
(5 > is an anomalous exponent depending on the initial data, and U2 decays to 
zero with normal rate 0{t~2'j. The idea of the proof is to combine the a priori 
estimates for the decay of global solutions with the renormalization group (RG) 
method for establishing the self-similarity of the solutions in the large time limit. 
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1 Introduction. 



In this paper, we study the initial value problem of the following thermal-diffusive com- 
bustion system: 

Ul,t = Ui^xx-Uiul (1.1) 



dU2xx + Uiul, xeR^, (1.2) 



with non-negative initial data {ui,U2)\t=o = {cli{x) , a2{x)) E {L^{R^) H L°°{R^)Y, of 
arbitrary size, where the positive constant d is the Lewis number. We are interested in 
the large time behavior of solutions of (|1.1| )- ([L^ ). 

The system (|1.1|)-([TT^) on a bounded domain is well-studied in the literature, see 



[Tp, IQ, [10 1, 1 11] and references therein. In case of homogeneous Dirichlet or Neumann 
boundary conditions, the large time behavior of solutions is that (mi,M2) converges uni- 
formly to a constant vector (ci, C2) such that ci ■ C2 = 0, see K. Masuda |1TT|| . 

More recently, the system ( |l.l|) -( pr2D on the line R^ has been proposed and investi- 
gated as a model for cubic autocatalytic chemical reactions of the type A + 2B 3B, 
with rate proportional to where Ui and U2 are concentrations of the reactant A, and 
the autocatalyst B. We refer to the interesting papers by J. Billingham and D. Needham 
@], 1^, for details. In 0] and the authors established the existence of traveling front 
solutions rigorously by shooting and phase plane methods; moreover, they studied the 
long time asymptotics of solutions by formal methods and numerical computations for a 
class of front initial data, i.e. data such that ai + 02 has nonvanishing limits as x — 00. 

Motivated by thermal-diffusive models with Arrhenius reactions, [|12], etc. Berlyand 
and Xin |^ considered system ( p..lD -( |]r2D for a class of small initial data in (L^nL°°(i?^))^ 
and showed that Ui{i = 1,2) are bounded from above and below by self-similar upper 
and lower solutions. The results of imply that Ui decays to zero in time with an 
algebraic rate faster than t^2^^^ for some 6 > 0, and U2 decays to zero like 0(t~2). 

In the present work, we prove the exact large time self-similar asymptotics with no 
restriction on the size of initial data as long as the data has sufficiently fast spatial decay. 
Our main result is the following. We consider the system ( p..l|) -( p!^ with initial data 
(01,02) E B X B, where B is the Banach space of continuous functions on R^ with the 
norm 

11/11 = sup |/(a;)|(l + |x|)^ with q > 1 fixed below. (1.3) 

xeR^ 

Let = 0(x) be the Gaussian: 

1 

</'(3^) = ^7|^e 4.. (1.4) 

Given A > 0, let ipA be the principal eigenfunction (ground state) of the differential 
operator: 
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on L'^{R^, dfx), with djj,{x) = e"* dx. The corresponding eigenvalue is denoted by Ea > 
(and Ea > for A > 0). We normahze ipA by / ?/'^(a;)(i/i(x) = 1. Our main result is: 



Theorem 1 (Global Large Time Self- Similarity). Consider initial data (01,02) G i3 x 
B,ai ^ OjOi > 0,i = 1,2. Let A = /^i ai{x) + a2{x)dx, the total mass of the system, 
which is conserved in time. Then system ( \1.]\ )- ^T7^ ) has a unique global classical solution 
{ui{x,t),U2{x,t)) G BxB for^t > 0. Moreover, there exists aq{A) such that, if q > q{A) 
in there is positive number B depending continuously on (01,02) such that: 

||t^+^^Mi(v^-,t)-i?^A(-)ILrO' (1-6) 

tjoo 

||t^W2(Vt-,t)-A0(-)|I^O. (1.7) 



Remark 1.1 All the results of Sections 2 and 3 hold for any g > 1. We need the decay 
at infinity of Oj to be fast enough only to obtain the exact decay rate in ( |1.6| ). For A 
large, Ea will be large, and the decay in ( |1.6| ) may be much faster than the diffusive 
one. Alternatively, we could consider data Oj G ^Bexp where ;Bexp is defined through the 
norm 

||/||,.p = sup|/(a;)|e^l^l 

X 

for some 7 > 0. Then, the conclusions of Theorem 1 hold for any A. 



Remark 1.2 The rate of convergence in ( |1.6| ) and ( |1.7| ) to zero is actually 0{t ^'), for 



some ?7 > 0, see ( |4.28| , [4.291) . The convergence in (|1.6| ) and (|1 . 7]) implies that 



ui{x,t) ~ — J — "T h h.o.t. and 

A X 

U2{x,t) ~ —(j){—=) + h.o.t., 
t2 y/t 

as t — > 00, where the leading terms are just the two parameter self-similar solutions to 
the reduced system. The anomalous exponent Ea occurs as a result of the interactions 
of nonlinearities of opposite signs. Furthermore, Ea can be computed or estimated as 
the ground state energy of operator Ca depending only on the Lewis number d and the 
total mass of the system. A nonperturbative upper bound is , while, for A small, 

Ea = ; O(A^), 

see P] for details. Actually, it is more natural physically to normalize A = 1, which 
amounts to putting a coupling constant A'^ in front of the reaction terms in ( |l.l| , p..2D , 
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and the anomalous exponent Ea depends then on the strength of that couphng constant. 



Remark 1.3 In order to understand the heuristics of ( |1.6|) , (|1.7|) , consider a more general 
problem: 

Ul,t = Ui^ccx - MlM™ (1-8) 
U2,t = du2,xx + UlU^ (1-9) 

for m > 1. For m > 2, as explained in 0, we can use the RG method of ||^ to prove 
that both Ui and U2 go diffusively to zero. For 1 < m < 2, one can use the maximum 
principle, as in Lemma 2.3 and equation ( ^.51 ) below, to bound from above Ui by ui, 
which is the solution of 

ui,t = Ui,.. - (1.10) 

Then using the Feynman-Kac formula, we get 

Mi(x,t) < exp(-C(t^-T)) (1.11) 

for |x| < 0{Vi). For |x| > 0{Vt), one gets a diffusive behaviour, depending on the 
rate of decay, as x ^ ±oo, of the initial data. Then, inserting the fast decay ( |1.11| ) of 
Ml in (|1.9| ), one shows that the effect of the nonlinear term in ( [1.9|) is small and that U2 
diffuses to zero. Clearly the borderline case m = 2 is the most delicate and the most 
interesting one. Instead of ( |1.11|) , one gets exp(— (9(logt)) which gives rise, after some 



analysis, to (KM) 



The rest of the paper is organized as follows. In section 2, we derive a priori estimates 
on the solutions of the system (|l.l|)-(pr2D based on the work of K. Masuda |1T[ for finite 



domains. Quite a few estimates are different here due to the unboundedness of R^. 
The a priori estimates imply the existence of global smooth solutions. In section 3, 
we derive decay estimates for the solutions using the maximum principle and a simple 
renormalization group (RG, see p|) idea to show that ui goes to zero like C(r 2"-^), for 
some 6 > 0. We use this information to prove that the nonlinearity is irrelevant (in the 
RG sense) in ( p.. 21 ) and that ||m2||oo < C>(t~2) as t — cx). In section 4, we use the results 
of sections 2 and 3, and the renormalization group method to prove the convergence to 
a self-similar solution and thus complete the proof of the main theorem. 



2 A priori estimates and global bounds. 

The goal of this section is to prove 

Proposition 1 The system ^1 . 1\ , has a unique classical solution satisfying 

WuiWhp < C{ai,a2)., 2 = 1,2, 1 < p < +oo, (2.1) 
where the constant C depends only on the initial data (01,02) G {L^{Ft-) fl L°°{R^)y. 
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Remark 2.1 Although some of the arguments below follow those of Masuda |Tl|, we 
provide them for completeness. Here and below, we use C to denote a generic constant 
that may vary from place to place. Moreover we write, as above, C(-) to indicate the 
only variables on which the constant may depend. 

First, we have the obvious 
Lemma 2.1 The solution (mi,M2) satisfies the L} estimates: 

11^1 + ■U2||Li(i?i) = + a2||Li(Ri) , ||M2||Li(i?i) > ||02||li(R1), 

\\ui\\l^(ri) < \\ai\\Li(Ri) , J uiuldxdr < +00. (2.2) 



Proof Integrating (|l.l|)-(pT^) over R^, assuming spatial decay at infinity, we get: 



Wih^it) = \\ai\\Li - [ \\uiul\\Li{r)dT, (2.3) 
Jo 

\\u2\\L^{t) = ||a2||Li + / \\uiul\\Li{T)dT. (2.4) 

Jo 



Combining (^)-(|23) gives (^. □ 



Lemma 2.2 The function gp{u2) = satisfies, for p > 2, 



0<9'M^)<{i^)9piu2)g;{u,)t^. (2.5) 



Proof Direct calculation. □ 

Using the classical parabolic maximum principle, we have. 
Lemma 2.3 

(1) <ui{x,t) <\\ai\\^, Vt>0; (2.6) 

(2) < U2{x, t) < U2{x,t), V t > 0, where M2 is a solution of: 

U.2,t = du.2,xx, lL2\t=o = a2{x)\ (2.7) 

(3) Mi(x, t) < ui(a;, t), Vt > 0, where ui solves: 

uit = ui^xx - ui ■ ml, Ui\t=o = ai{x). (2.8) 

Remark 2.2 Using this Lemma, one immediately proves ( |2.1| ) for ui, since, by ( f2.8D , Ui 
is a fortiori bounded by the solution of the heat equation with the same initial data. 

Lemma 2.4 The solutions (mi,M2) of (|rT])-(|T]^ satisfy the If bounds: 

\\ui\\ LP ^ C {0,1, a2,p) < +00, i = 1,2, 1 < p < +00, and p = integer. (2.9) 
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where C{ai, a2,p) is a constant depending only on the initial data and p. 

Proof Due to Remark 2.2, we have only to prove the bounds for U2- We use standard 
local existence of classical solutions for parabolic equations, and, therefore, we freely 
integrate by parts below. Our goal will be to prove bounds uniform in time. We shall 
show that 



/ U2dx + / / {ul U2 + ul U2 + UiU2^'^)dxdT 
JR^ Jo JR^ 

< C(ai,a2,p)(l+ I I uiul^^dxdr) (2.10) 
Jo Jr^ 



'0 JR^ 

for all j9 > 2 ( p integer). Besides, we shall show, for p = 1 

/ / (uiul + ul U2)dxdT < C(ai,a2)(l+ / / uiuldxdr) (2.11) 
Jo Jr^ ' Jo Jr^ 



Using ( p.2|) to bound /q /^i Uiuldxdr, and using induction in p, we get that all the terms 
on the left hand side of ( |2.10|) are bounded, for all p > 1, p integer {u2 G by (pl^)). 
In particular, this implies the claims of the Lemma. 

To prove ( |2.10| , |2.11| ), we multiply (|L2|) by (7^(1x2), we integrate over x Ri, and 
we get, after integrating by parts: 



U2)dx = j ^ gp{a2)dx ~ d J J ^ u\^^g'p{u2)dxdT 



+ [ [ uiu\g'iu2)dxdT, (2.12) 
JO Jr^ ' 



here gp{u2) = u^, p> 1. 
Consider the identity (p > 1), 

^ J^^iuj + Ui)gp{u2)dx = J^^{1 + 2ui){ui^xx - Uiul)gp{u2)dx 

+ {ui + u\)g' {U2){du2.xx + uiul)dx 
Jr^ ^ 

{I + 2ui)ui^xU2,xg'p{u2)dx - 2 / ul ^gp{u2)dx - d I {I + 2ui)ui^xU2,xg'p{u2)dx 
_Ri ^ Jr^ ' Jr^ ^ 

d I [ui + u^g'p{u2)u\,j.dx — I {1 + 2ui)uiu\gp{u2)dx + I {ui + u^g'p{u2)uiu\dx 
I + 11 + in + IV + V + VI. (2.13) 



We estimate for p > 2, using (p^), 

/ + // + ///< (1 + 2||ai||oo)(l + rf) / \ui,xU2,x\gL{u2)dx - 2 u\^gp{u2)dx 

Jr^ ^ Jr^ 

< (1 + 2||ai||oo)(l + rf) / \ui.xU2,x\i9piu2)g''iu2)i -))^dx 

Jr^ ' ^ p — 1 
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-2 / , ul.^9p{u2)dx 
< l(l + 2||ai|U)(l+rf)(-3Y)^(e I ulMn2)dx + e-' f ulj;{u2)d 



-2 / ul^gp{u2)dx (2.14) 

Picking 

' {l + 2\\a4^){l + dy p ' 

in (2.14), we continue: 



I + 11 + III 



< (^)e-i(l + 2||ai|U)'(l + rf)' / ulX{u2)dx- f ul,g,{u2)dx (2.15) 

In addition, we have: 

IV < 0, (2.16) 

V < — Uiulgp{u2)dx, (2-17) 

VI < (||ai||oo + ||ai||^) / g'p{u2)uiuldx. (2.18) 



Integrating (|2.13|) from zero to t yields: 



< ( / {ui + uj)gp{u2)dx){t) < (ai + aj)gp{a2)dx 



+C{p,ai) f f ul g''{u2)dxdT 
Jo Jr^ ^ 

+ (Ikilloo + llflillL) / / g'{u2)uiuldxdT 



Jo JR^ ^ 



J J uiulgp{u2)dxdT — J ^ul^^gp{u2)dxdT. (2-19) 



Combining (g]T2|) and ( p39D gives (|2J^ ) for p > 2. 



For p = 1, we proceed from (|2.13|) as follows: 

— / (Ui + u'^)u2dx = — [ (1 + 2ui)ui xU2 xdx — 2 [ ^M2dx 
dtJR^ Jr^ ' ' ' Jri ^'^ 

—d / {1 + 2ui)ui xU2 xdx — / {1 + 2ui)uiuldx 

JR^ ' ' JR^ 

+ (ui + u\)uiu\dx 

JR^ 

< (1 + 2||ai||oo) J^^ \ui,xU2,x\dx - 2^^ ul^^U2dx 
+d(l + 2||ai||oo) / \ui^xU2,x\dx - / uiuldx 

JR^ JR^ 



+ ( Ctl oo + 



killoo) / Uiu\dx 
Jr^ ' Jr^ 



+ (||ai||oo + / uiujdx, 

R^ 



(2.20) 



for any e > 0. Now, integrate (|2.2CI|) from to t, to get (for e small enough): 



/ (2uiul + ul U2)dxdT < 2d / / u\ dxdr 
' Jo JR^ 



+ C(ai,a2)(l+/ / {ul^ + UiuDdxdr). 
Jo Jr^ 

Now, use ( |2.12| ) with p = 2 to bound 

2d / ^2 x^^d'^ — / ct^dx + / / uiuldxdr. 
Jo Jr^ ' Jr^ Jo Jr^ 

Finally, observe that, multiplying ( |1.1| ) by Ui, and integrating by parts, we get: 

— / u'idx = — [ M? ^dx — I u^u^dx 

dtJRi jRi ^'"^ Jr^ 



IR^ JR^ 

from which we immediately obtain: 

rt 

lo JR" 



(2.21) 



(2.22) 



(2.23) 



I I u\ dxdr < / afdx 
Jo Jr^ ' Jr^ 



uniformly in t. Combining (|2.21| , |2.22| , p.24| ), we get ( |2.11|) . This completes the 
the lemma. 



(2.24) 

proof of 
□ 



Remark 2.3 In Masuda [|TT|, gp{u2) = (1 + U2Y1 and lemma 2.4 is proved for fractional 
p's by starting the induction from ( |2.12| ), p G (0,1) instead of combining ( |2.20D and 
( p.l9| ). However, for unbounded domains like such argument fails since gp{u2) ^ 
L^{R^). Lemma 2.4 still holds for fractional powers p, and we will show that with the 
help of L°° bounds which we establish below using ( p.lO| ). 



Lemma 2.5 The solutions (mi,M2) of ( |l.lD -( pr^ obey the estimates 

\\ui,xh + \W2,xh < C{ai,a2,ai^x,a2,x) (2.25) 
where C{ai,a2,ai^x,(^2,x) is a positive constant depending on ||aj||i, HctiUcxD, llc^j.xlh- 

Proof Multiplying ( |1 . 1|) by Ui^xx, (|1-2| ) by U2,xx, and integrating over gives: 

1 d 



2dt 



hi,x\\l 



/ \uixx\'^dx— / Uiuiuixxdx 
Jr^ ' Jr^ ' 
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or 

1 II ||2_ 

Similarly, we have: 
Id 2 



/ \uixx\'^dx— / MoU? „(ia; — 2 / U2UiUixU2xdx (2.26) 



-(i / \u2xx\'^dx + 2 / I 'u\uixU2xdx. 

JR^ ' ifli 'J ' ' 

Adding (^1261) and (|2:2^ ) gives: 



(2.27) 



I'^/ii ||2 II ||2 
2^(I|WmII2+ \\U2,x\\2 



— I UnU, ^dx 
Jr^ ' 7i?i ^'"^ 

2 y UiU2Ui^xU2,xdx + 2 y MlM2W2,^(ix + y M2Ml,a;M2,a;C^a;. (2.28) 



Now, integrate (|2.28|) from to t, and bound the resulting terms on the right hand side 
as follows. The first three terms are negative, and the last three terms can be bounded, 
using 11^1 lloo < Ikilloo (see (p.6|)) and the Cauchy-Schwarz' inequality: 



/ / U1U2U1 xU2 xdxdrl < ||ai||oo(/ / u\ uidxdr / / u\ dxdrY^'^ (2.29) 
JO JiJi ' ' Jo JR^ ' Jo Jr^ 



/ / UiU2ul dxdrl < ||ai||oo / / U2u\ dxdr, 
Jo Jr^ ' Jo Jr^ 



(2.30) 



and 



/ / ului xU2 xdxdrl < (/ / u\u\^dxdT / / u^u^^dxdrY^'^ . (2-31) 
Jo Jr^ ' ' Jo Jr^ ' Jo Jr^ 



Now, the terms on the right hand side of ( p.29| , p.30| , p.31|) are uniformly bounded in t, 
because all the terms in the left hand side of ( |2.1U| ) are bounded, for p >2. So, we get 



I^*i,2^ll2(^) + ll^*2,x||2 W < C{ai, 02, ai^x, a2,x), V t > 0. 



(2.32) 
□ 



By local existence of classical solutions in C([0,T*); {L\R^)f) n C((0,T*); {H\R^)f) 
of system ( |1.1|) -( |L2|) , and parabolic regularity: 3ti > 0, such that 

\\ui,x\\l + ll^^2,x||2(ii) < C{ti){\\ai\\l + ||a2||2), h < T* 
If we replace by ti in the proof of Lemma 2.5, then (|2.25| ) in fact implies 
Corollary 2.1 

\\ui,xh + \\u2,xh<C{ai,a2) (2.33) 
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where the constant C depends only on ||aj||i and \\ai\\ooi'i = 1,2. 
Corollary 2.2 By Sobolev imbedding, ( ^.!^ ) forp = 2 and \2. 3S[) imply: 

IIM2II00 < C(ai,a2). 
Thus Proposition 1 follows from: 



(2.34) 



Corollary 2.3 Combining the estimates in Lemma 2.4, ( ^-dj ), ( 2.34 ) o-nd standard local 
existence of classical solutions, we have shown that global smooth and bounded solutions 
exist for system f^-^^ m (L^R^) n L'^{R^)f. 



3 Decay estimates 

To exhibit the decay of the solutions of (fL.l|)-( pr2D , let us introduce the scaled solutions 

■Ui{x,t) = Viui{Vix,t) 

for i = 1, 2. From now on, we consider nonnegative initial data (ai(x), 02 (x)) ^ B x B 
{L^{R^) n L^{R^)Y for g > 1. The purpose of this section is to prove 

Proposition 2 The solution {ui,U2) of constructed in Proposition 1 satisfies the 
bounds 

\\ui\\{t) < C(||ai||,||a2||)(l + t)-' (3.1) 

\\u2m < C(||ai||,||a2||) (3.2) 



where S = 8{\\a2\\L^) > ^ and || • || is, for all t > 0, the norm (l.S). 



Remark 3.1 Note that, in particular, (|3.1| ) and ( |3.2|) imply 

Ikilloo(t) <C^(||ai||,||a2||)(l+t)-^-^ 

||M2||oo(t) <C(||ai||,||a2||)(l + t)-^ (3.3) 

Remark 3.2 The bound ( ^^TD was essentially derived in [H, using ( p.8|) and ( p^) , but 
with a different norm. 



Proof Using (|2.^ ), (|2.7| ), it is enough to prove ( |3.1|) with ui replaced by Ui. Also, since 



U2 solves the heat equation, we have: 



M2(x,t) = ^0(^) + /l.O.t. (3.4) 
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as t — s> oo, with given by ( |1.4| ) and A = J^i a2{x)dx. 

The higher order terms are easy to control, so we can consider, instead of ( p.8| 



ui,t = - (3.5) 
By a simple change of variables, ^ = r = logt, we get, for t > 1, 

Vm{x,t) = (e"^^-^Mi(-,l))(0 (3.6) 



with Ca given by ( |1.5| ). Then, ( |3.1| ) follows from the first bound on the semigroup e"'^^-'^ 
given in Lemma 4.1.b. 
To prove (|3.2|), we need 



Lemma ^.1 (Sharp Decay of U2 in L°° norm). There exists a constant C depending only 
on \\ai\\, i = 1,2 such that 

||n2||oo(t) < —^-T, Vt > 0. (3.7) 

(l+t)2 



Consider then equation ( |1.2|) : U2,t = du2,xx + {uiU2)u2- By the first inequality in ( p.3|) 
(which follows from ( p.l|) ), and Lemma 3.1, U2 < ^2, where U2 solves: 

U2,t = dU2,xx + C{1 + ty^'^U2 
U2\t=0 = «2- (3.8) 

But ( p. 81) is a linear heat equation, from which (|3.2| ) follows easily since {l+t)~^~^dt < 
oo. □ 

We are left with the 

Proof of Lemma 3.1 Write equation (|1.2|) in integral form: 

U2{t,x) = / - e i'ita2{x - y)dy + 

\ (47r(is)^2 exp( —){uiul){x — y,t — s)dyds. (3.9) 

JR^ Ms 

Taking L°° norm in ( p.9|) yields 

||M2||oo(t) < ^(11 ^2 II) rs-^||«iM^||Li(t-s)rfs (3.10) 

(l+t)2 ^0 

Now, use ||'Ui'U2lUi — II'^^iIIooII'wIIIlij ||wi||oo(^) < C'dkill? Ik2||)(l + 1)^^^"^^^ (which is 
the first inequality in (|3.3|)), 



\ul\\L^{t) < \\u2\\oc{t)\\u2\\L^{t) < C{\\ai\\, llaall), (3.11) 
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(which follows from Proposition 1) and 

r + {t- s))-^/^-^ds < C(l + t)-^ (3.12) 
Jo 

to get 

h2||oo(t) <C(||ai||,||a2||)(l + t)-' (3.13) 



(for 6 < 1/2, or else, (p.7| ) is already proven). Now, use (|3.13|) to improve ( p.llD into 
||Mil|Li(t) < C'dlaill, ||a2||)(l + t)~^ which, inserted into (|3TT0| ) yields 

h2||oo(t)<C(||ai||,||a2||)(l + t)^'^ 

and we can iterate up to 1/2, which proves (p.7|). □ 

4 Self-similarity 

In this section we apply the Renormalization Group method to improve Proposition 2 
and finish the proof of the Theorem. We prove: 

Proposition 3 Under the assumptions of Theorem 1, there exists e > such that, if 

||ai||||a2|| < e, (4.1) 

the claims of Theorem 1 hold. 

Proof of Theorem 1 By Proposition 2, we can find a T < oo such that the functions 

UiT{x,t) = \/Tui{\/Tx^Tt) satisfy 

lkiT(-,l)||lk2T(-,l)ll <e, (4.2) 

where T depends on the initial data (ai, 02). Moreover, UiT solve the equations ( |1 . 1|) and 
( p.. 21) , and thus, by Proposition 3 and (|4.2| ), UiT and hence Uj, will have the asymptotics 
claimed in the Theorem. □ 

We will now set up an inductive scheme for the proof of Proposition 3. We define, 
for L > 1, 

u^^\x, t) = U'uiiVx, L'-^t), t e [1, L"]. (4.3) 
Then ■Uj-"'' satisfy the equations ( |1 . 1| ) with initial data 

at\x) = U'ui{L''x,L^''). (4.4) 
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We will study a!f^ inductively in n i.e. we will consider the RG map (ai, 02) — > (a'^, a'g) 
defined in B ^ B where 

a'iix) = Lui{Lx, L^) (4.5) 



and {ui,U2) solve (|l.l|) -(pr^ ) with initial data (01,02). We first prove a Lemma for the 



linearization of this map when ( |l.l| )-( pr2D is linearized around the expected asymptotics 



§ and (jlj) 



Hence we consider the equations 



vit = vixx - AH~^(j)'^{--^)vi 

V2t = dV2xx (4.6) 

for t G [1, L^] and Vi{x, 1) = aj(x). By the change of variables ^ = ^^'^ = log^ gets 

Lwi(Lx,L2) = (L-2^-^ai)(x) 

Lv2{Lx,L') = (L-2%)(a;) (4.7) 

where 

I d 1 , , , 

C = -d— - -X- 4.8 

dx^ 2 dx 2 ^ ^ 

and Ca is given by (|1.5|) . Recall also that ipA is the principal eigenvalue (ground state) 
of Ca- We collect some properties of Ca and C in 

Lemma 4.1 Ca and L have the following properties 
a) Ca^a = EaiPa. Ea>Q z/ a > 0. 

h) Let f E B. There exist 6 > and tq < 00 such that, for r > tq, 

||e--^-/|| <e-^i/||. 

Moreover, there exists q{A) such that, if f E B, with {ipA, f) = ({■,■) being the 

^2 

scalar product inH = L'^(R, djj), dji = e~dx), and with q > q{A) in ( \T7^ ), 



X 



c) Let Pa be the orthogonal projection in Ti on ipA- The quantities \Ea — Ea'\, |1 — 
{ipA,'ipA')\, 1 1 -Pa — Pa'W (operator norm in B) are bounded by C{K)\A — A'\ and 
\\Pa\\ < C{K), for 0<A,A'<K. 

d) = 0. 

e) Let f E B, J fdx = 0. Then, there 6 > and tq < 00 such that, for r > tq, 

||e--^/|| <e--^' 
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The proof of Lemma 4.1 is based on [0 (see also [T^), and will be given in the 
Appendix, where we also show that the scalar product {ipA, f) is well-defined for f E B. 

Returning to the proof of Proposition 3, we write the RG map, defined by ( [4.5|) , as 

a[ = L~^^^ai + Lni{Lx,L'^) (4.9) 
= L'^^a2 + Ln2{Lx,L^) (4.10) 



where 



ni{x,t) = - I ds f dyGAit,s,x,y){uiul{y,s) -UiA^s ^(f)'^{^)) (4.11) 



n2{x,t) = J ds J dyG{t — s, X — y)uiul{y, s) (4-12) 



and Ga is the fundamental solution of the vi equation in ( |4.6| ) and G is the kernel of 
gd(t-s)A^ where we write A for Denote by sl the scaling (sl/)(x) = Lf(Lx) and by 
GA{t,s) the operator corresponding the kernel GA(t, s,x,y). Then we have 

Lemma 4.2 

a) \\sl\\ < L 

b) \\GAit, s)\\ < e^(*-^); ||e'^(*-^)^|| < e'^^*"^), for c < 00. 
Proof a) follows from 

supL \f{Lx)\{l + \x\y <L\\f\\ 

X 

and b) from 

< GA{t, s, X, y) < Go{t - s,x -y) 

(which itself follows from the Feynman-Kac formula and ^4^0^ > 0), and the explicit 
Gaussian kernel of Gq. The kernel of e'^^*"'')^ is similar. □ 

Let us now specify the A in ( [4.9|) and ([4.11|) (which is not the same as the one in 
Theorem 1). We write 

ai(x) = BipAix) + bi{x) (4.13) 
02 (x) = Acpix) + b2ix) (4.14) 

with 

B = {i'A,ai), 



A = y a2dx. 
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Remembering that ipA is normalized by (ipAji'A) = 1, and by J (p{x)dx = 1, we see 
that 

j b2dx = 0. 

Write {a'i,a2) similarity, with primes. The main estimate then is 

Lemma 4.3 Given L > Lq = e^"^", with tq as in Lemma 4.1, there is an eo{L) > such 
that if ||ai|| ||a2|| < e < eo{L), we have 

a) \A' -A\ < C(L)e||a2|| 

b) ll^'all < L-^^\\b2\\+C{L)e\\a2\\ 

c) \B' - L-^^-^B\ < C(L)e[e(l + ||a2||) + II62II] 

d) \\b[\\ <L-2(^^+'^)||6i||+C(L)e[e(l + ||a2||) + ||62||] 
where C{L) is an L-dependent constant. 

Proof We solve first U2 from the equation 

U2ix, t) = e'^(*-^)^a2 + n2ix, t) = U20 + n2, (4.15) 
with n2 given by ( [4.12| ), by the contraction mapping principle. Consider the ball 

Br = {u2 : \\\u2\\\ = sup ||m2(-, ^)|| < -R||a2||}- 

te[i,L2] 



For U2 G Bf( we bound n2 by using (|2.6|) , i.e. 

||«i(-,s)||oo < lle'^^^-^^Vlloo < llailloo <C||ai|| 
and Lemma 4.2.b, to get 

IIN2III < R^C{L)\\ai\\ \\a2f < R^C{L)e\\a2\\. (4.16) 
Since by Lemma 4.2.b again, 

|||^t2o||| < C'(-^^)||a2||, 

we see that the right hand side of ( [4.15| ) maps Br into itself if i? = R{L) is large enough 
and e < e(L) is small enough. It is easy to see that the right hand side of ([4.15|) is a 
contraction in Br, so that we get a solution in Br. By Lemma 4. 2. a then, 

\\Ln2{L;L')\\<C{L)e\\a2\\ (4.17) 

and since 

A' = A + J Ln2{Lx,L'^)dx 

h'2 = L-^% + Ln2{Lx,L^) + {A-A')(P, 



15 



a) and b) follow from ( [4.17| ), and Lemma 4. I.e. 



For a'l^ consider rii in ([4.11|) , and write 



w{y,s) = ul{y,s) - A'^s 

= My, s) + A.-i/20(^))((e'^(-i)^62)(x) + n,{x, s)), (4.18) 

using (CT) , (CT ) and e^("-i)^0(y) = s-^{^). Thus, by (CT) and Lemma 4.2.b, 

llhlll <C(L)(||a2||+A)(||62||+e||a2||) (4.19) 

and, since A < C||a2||, sltli = Lni(L-,L^) is bounded by 

\\sLni\\<C{L)e{\M+e\\a2\\). (4.20) 

(use ||ui||oo < C'llaill and ||ai|| ||a2|| < e). 
Since, from 



B' = (V'A', a'l) = {i^A',L ^^^ai) + {ipA', slUi) 

and, from ( |4.13|) 



(we used Pa&i = 0), we get, using Lemma 4.1.c and ([4.20|) , 

\B'-L-'^-B\ < C\A-A'\{B+\\h\\) + CiL)e{\\b2\\+e\\a2\\) 

where the constant C in C\A — A'\ is independent of A, A', because A here is uniformly 
bounded (by Lemma 2.1). Then, using part a) above, we get 

\B'-L-'^-B\ < C(L)e[(i?+||6i||)||a2|| + ||62||+e||a2||] 

< C(L)e[e(l + ||a2||) + ||62||)] (4.21) 

(since i? + < C||ai|| and ||ai|| ||a2|| < e), i.e. we prove c). 
Finally, for b[, write (use {J^, (|^) 

b[ = (1 - PA')a[ = BL-^'^^iPA - Pa')^a + L-^^^b^ + (Pa - PA')L-^''^b^ 

+ (1 - PA')sLni 

(using again Pa&i = 0). Now, Lemma 4.1.b, c and ( [4.2UD imply 

||6;||<L-^(^-+^)||6i||+C(L)e[e(l + ||a2||) + ||62||] 

which is d). □ 

For later purposes, we derive a lower bound for B. Recalling the definition (4.18), 
write 

uu = (A - aW(^))wi - {w{x,t) + {A' - A'')r'<f)\^))m. 
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Using the Feynman-Kac formula, we deduce 

a\{x) > (L-2^^'ai)(x)e-^(^)(lll-'lll+l^'-^"l) 
and thus a lower bound 



Proof of Proposition 3 We decompose as in ( |4.13| ), ( |4.14| ) and derive bounds for 
An, Bn and 6" using Lemma 4.3. Set 



nE„ 



n-l 
m=0 



Note that > A, so that Ea„ > -E'a > 0. Let r] < min(5, E'^)- Then there exists a 
constant C{L) (depending possibly on L but not on n) such that 



< A„ 


< 


C{L)\\a4 


Ir2 II 


< 


C(L)L-2"''||a2|| 




< 


C(L)L-2"^"||ai|| 


ll^ill 


< 


C(L)L-2"(^"+'')(||ai|| +e||a2 


en 




IKIIII«2ll < C(L)L-2"^"e. 



(4.23) 

The bounds ( |4.23| ) hold by definition for n = 0, and the induction follows from Lemma 
4.3: the bound on follows from the first four bounds in ([4.23|) , and it can, in turn, 
be inserted in Lemma 4.3 to iterate those bounds. For i?„, we iterate Bn < C{L){1 — 
j^-nrj^ j^-2nE„ ||^^ || (^-^j^jcj^ implies (|4.23| )), in order to control the right hand side in Lemma 
4.3 c. Furthermore, the bound on e„ and Lemma 4.3 a imply that 

\An+i - A^\ < C(L)L-2"^"e||a2|| (4.24) 

and thus An ^ A*, for some A*] moreover, 

A* = J (ai + a2)dx, 

because /(ai + a2)dx is conserved (by Lemma 2.1) and / a^dx 0, by ( [4.13| ), ( [4.23| ). 
Since Ea is continuous in A, by Lemma 4.1 c. 



TP . TP* 77' Z? ^ Z?* 

From Lemma 4.3.C and ( ^4.23| ), we get that 

\Bn+i - < C(L)eL-2"(^"+'')[ 

and by ( [4.25| ), there exists a B* such that 

S^L^"-^" ^ B* 



P2 



(4.25) 



(4.26) 



(4.27) 
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By and (^), (Q), ( pl , 



so, B* > 0. Equations ( PB[ ), ( ^IQ and ( ^71) may be rewritten, using (|^, 
(ill), as: 



\\Viu2iVt-,t)-A*(P\\ < Ct-"! 



0-2 



(4.28) 



(4.29) 



for times t = L^", L > Lq. For t G [-Z^^", L^""*"^] we use similar estimates for the rii in 
(|]T3) and (^4.121) , and, dropping the *, we get (|L6|, [l^)- □ 
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Appendix 1: Proof of Lemma 4.1. 

First, observe that Ca, acting on its domain in L'^ {R^ , dfi) , is conjugated to a per- 
turbation of the Hamiltonian of the harmonic oscillator: 

2 2 (7^ T'2 1 

c-Cac-- =Ha = + - i + ^V(a:), (A.l) 

acting in L'^{R^,dx). Hence, Ca has a compact resolvent, a pure point spectrum and, 

using the Feynman-Kac formula and the Perron- Frobenius theorem , a non-degenerate 
lowest eigenvalue. The same conclusions hold for C 
To prove a), let us differentiate 

CAi^A = EAi^A {A.2) 

with respect to A^. We get: 

+ Ca^^'a = E'a^^a + Ea^P'a- {A.3) 

Now, we take the scalar product of (A. 3) with ipA, and use {ipA,'ipA) = 1 (which implies 
(V'a,^a) = 0), to get 

E'a = {^a,<P'^a). (A4) 
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0, we have ipo 



.2 



/47r 



and En = 0. Therefore, 



Since > 0, E'^ > 0, and, for A 
EA>OhT A> 0. 

To prove b), we discuss only the second claim, since the first one is similar but easier 
(and holds for any g > 1). Observe that, since Ea is non-degenerate and / is orthogonal 
to ^pA, the bound would be trivial if we took the norms in Ti. But functions in Ti. have 
essentially a Gaussian decay at infinity, while those in B have a polynomial decay. To go 
from a contraction in 7i to a contraction in B, we use an idea of [0. First notice that, 
since y4^0^ > 0, the Feynman-Kac formula gives 



-tCa 



and e '^^°{x,y) is explicitly given by Mehler's formula [13 
(e-^^°)(x,?/) = (47r(l - e-^))-5 exp ( 



X — e 



4fl 



(A.5) 



(A6) 



Hence, if a function v satisfies 



for some constant C, we have 



\V[X] 



<C(l + |x|)-^ 



|(e"^^"t;)(x)| < C'e^{l + \x\e^y 



{A.7) 



{A.8) 



for |a;| > 2y/qf and another constant C Hence, the operator 6^"^^° contracts, for |a;| 
and r large, any function that decays as in (A. 7) with g > 1. By (A. 5), we see that Ca 
behaves similarly. So, to prove b), we shall use the contraction in Ti, for x small and 
(A. 8) for X large. However, we need here q large, depending on Ea, hence on A. For the 
other bounds in Lemma 4.1, any g > 1 suffices. 

Besides, let 0„ be the nth Hermite function which is an eigenvector of Hq in (A.l) 

(they are of the form Pn{x)e s , where P„ is a polynomial of degree n). One can show 
that, for any C > 0, for some even n = n{A) and for any |x| large enough. 



iHA-EA)iC(Pn-e'T^lJA)>0. 



Indeed, {Ha — EA)e 



8 ipA 



by (A.l, A.2) and Ha^ 



> 



> (because 



> for 

n even and |x| large). Using the maximum principle, the inequality fg ~ | ~ Ea > 0, 

for \x\ large, and the fact that there exists a large |x| so that C0„ — e s > 0, for a 
sufficiently big C, one concludes that ipA is bounded by: 



< i^Aix) < C{A){1 + |x|")e- 



4 



{A.9) 



for some n = n{A), which implies that the scalar product {'ipA, f) for / G i3 is well- 
defined, if g = q{A) is large enough. 
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To prove b), it is convenient to introduce the characteristic functions 

Xs = x(kl < p) 
Xe = x(\x\ > p) 

where p will be chosen suitably below. The properties of Ca that we need are summarized 
in the following 

Lemma A.l There exist constants C < oo, c > 0, such that 

i) For g eB, 

||e"^^^fif|| < Ce^ll^ll. (AlO) 

ii) For g e L'^{R^,dn), 

||e-^-5ll < C\\g\\2, (All) 
where || • II2 is the norm in L^{R^, d/i). 

iii) For g such that Xsg £ dfj.), 

\\xee-'''^Xs9\\<e-^\\Xs9\\2, (A12) 

for p large enough. 

iv) For g e B, 

\\Xee-''''^g\\<e-'''^qg\\, (A13) 
for p large enough, and g > 1. 

Now, take f ^ B, {ipA,f) = 0, ||/|| = 1. We set r„ = np, and, using the Lemma, 
we prove inductively that there exists a. 6 > such that f (r„) = e''^"'-'^f satisfies, for p 
large, 

\\Xsv{Tn)\\2 + \\Xsv{Tn)\\ < e^e'^'^, {A.U) 

and 

\\xMrn)\\ < e-f'^ (A15) 

with P = [Ea + S)p. Part b) of Lemma 4.1. follows from (A. 14), (A. 15) by taking a 
smaller 5, in order to bound the constants, for r > Tq, with tq large (for times not of 
the form r = np, use (A. 10)). The bounds (A. 14), (A. 15) hold for n — 0, for p large 
enough, using ||/|| = 1 and the obvious inequality 

||Xs/||2<e^||/||. (A16) 

So, let us assume (A. 14), (A. 15) for some n>0 and prove it for n+1. Let v — i'(t„) 
and write 

V = XsV + XiV = Vs+ Vi. 
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For all n, (f (r„), ■j/'yi) = 0, so that = \ivi,ipA)\ < C{A)p ^'^ ^\ where we use 

(A. 9) to derive the last inequality. Then, we get 



4' 



where E'j^ is the second lowest eigenvalue of Ca and, in the second inequality, we choose 
S small, q > 6p and p large. For the third inequality, we used (A. 14) and p large, so that 
^ i_ Combining (A. 11) and (A. 17), we get 

again, for e~^P small enough. Finally, from (A. 10) and (A. 15), we have 

\\e-f^^Vi\\ < Ce^e"^" {A.19) 
and, from this and (A. 16), we get 

\\Xse"^^^ve\\2 < Ce^eh-^"". (A20) 



Combining (A.17)-(A.20), one gets (A. 14), with n replaced by + 1 for p large 
enough. On the other hand, (A. 15), with n replaced by n + 1, follows immediately 
from (A. 14), (A. 15) and (A. 12), (A. 13), taking cq > E + 6. We choose the constants as 
follows: take S small and p large and q > 5p > 

Turning to c), we observe that (A. 4) and ( |1.4|) imply that E'^ < (And)^^. Next, (A. 3) 
imphes that 

i;'^ = (Ca - Ea)-\E'^ - <\>^)i,A {A.2\) 

where 

{La - EaY^ = / e-"(^^-^^)dr (A.22) 
Jo 

is a bounded operator on the subspace {/ G B\{ipA, f) = 0}, because of b) above. Also, 
(A. 4) means that [ipA, {E'a — (f^)'4' a) = 0. Hence, the norm of the right hand side of 
(A. 21) is bounded, and we have 

Ua-^^'aW < C{A,A')\A'' - A'^\ < CiK)\A~A'\ (A23) 

for A, A' < K. Since (A. 9) shows that Pa is well-defined and bounded in B, point c) 
is proven. Point d) is an explicit computation (£0 = 0), and the proof of e) is similar 
to the one of point b), since / fdx = means that (/, 0) = where (■, ■) is the scalar 

2 

product in L^(_R^, e^dx) and 0, given by ( |1.4| ), is the principal eigenvalue of £. □ 

We are left with the 

Proof of Lemma Al. Part (i) follows immediately from (A. 5) and (A. 6). 
For (ii), we use the Cauchy-Schwarz inequality applied to 

{e-^^g){x)= / e-^^{x,y)e-i\^\{y)dy 
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and the bound 

/ f 2 \ 1/2 

which follows from (A. 5) and (A. 6). 

For (in) proceed as in (ii) by using Cauchy-Schwarz' inequality, but replace (A. 24) 

by 

sup(l + |a;|)^(/ \e-''''^{x,y)W{\y\ < p)e-"^"dyy < e'^^ (A25) 

\x\>p \J / 

which again follows from (A. 5) and (A. 6) (we can replace | in (A. 25) by | — e for any 
e > 0, if p is large enough). 

Finally, (iv) follows from (A. 5) and (A. 8). (Since it is enough to have q > 5pwe have 
\x\ > p > '2..^/qp for S small and we can use (A. 8)). □ 
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